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Abstract. We report new experimental results obtained on the insulating spin glass CdCr1.7Ino.3S4. Our experi- 
mental setup allows a quantitative comparison between the thermo-remanent magnetisation and the autocorrelation 
of spontaneous fluctuations of magnetisation, yielding a complete determination of the fluctuation-dissipation relation. 
The dynamics can be studied both in the quasi-equilibrium regime, where the fluctuation-dissipation theorem holds, and 
in the deeply ageing regime. The limit of separation of time-scales, as used in analytical calculations, can be approached 
by use of a scaling procedure. 
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1 Introduction 

Despite their large diversity, glassy systems have many dynam- 
ical properties in common. In particular, a similar ageing be- 
haviour can be observed in polymers, gelatins, or spin glasses 
11121 . Stationarity cannot be reached in these systems in ex- 
perimental, or even in geological times: they always remain 
out-of-equilibrium, even when not submitted to any external 
perturbation. 

During a long period, the theoretical activity was concen- 
trated on the study of the statics of glassy models. With the 
nineties, began the time of theoretical dynamical studies, first 
by numerical simulations, and then by analytical results on 
specific mean-field models 13 1415 1 . From these studies, new 
concepts appeared, generalising the well known Fluctuation- 
Dissipation Theorem (FDT), which holds for equilibrated sys- 
tems 16171 . 

In equilibrated systems with time translational invariance 
(TTI), FDT can be used to measure the temperature in an abso- 
lute way: 

d tw C{t-t w ) 



k B T = 



R(t - U 



(1) 



In this relation C(t — t w ) is the autocorrelation function of an 
observable (for instance the magnetisation M (t)) between two 
times, t and t w , and R(t — t w ) the response function associated 
with a pulse of the conjugate field at time t w , h(t) = S(t — t w ), 
These quantities are two times quantities, but, as the system is 
TTI, they depend only on the time difference t — t w . 
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Spin glasses never reach equilibrium, and the time autocor- 
relation and the response function can not be reduced to one- 
time quantities. Therefore, the temperature cannot be defined 
on the basis of usual concepts. Nevertheless, it has been shown 
that, in specific models with low rate of entropy production, 
and using a generalisation of the FDT relation, a quantity that 
behaves like a temperature could be defined [5 1, the "effective 
temperature". The effective temperature for one given value of 
C(t w ,t) = C can be defined as: 



knT f 
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lim 

t w — >oo 
C(t w ,t)=C 
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(2) 



The only difference between relations[2and|2]concerns the do- 
mains of validity. The generalised fluctuation-dissipation rela- 
tion is valid for stationary systems (simply, C and R depend 
only on t — t w and T e f / = T), and it is also valid for every 
systems in the limit of small rate of entropy production. Glassy 
systems, in the limit of long waiting time are such systems. 
Some experiments have been set up to measure this effective 
temperature using frequency measurements in glassy systems 

IMS). 

To understand the meaning of the time limit in equation^ it 
is helpful to refer to the so-called "Weak Ergodicity Breaking" 
(WEB) concept 11 II . WEB was introduced first in the study of 
the dynamics of a random trap model very similar to the Ran- 
dom Energy Model (REM) 1121 . According to WEB scenario, 
two different contributions can be identified in the dynamics: 
a stationary one, corresponding to usual equilibrium dynamics 
in a metastable state (and then not relevant for ageing studies), 
and a second one, describing the long term evolution between 
many metastable states, which features the ageing properties. 
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This approach agrees well with an experimental fact: in 
glassy systems, the relaxation function can be decomposed in 
two distinct contributions 1131 : 

- The first one is independent of the age of the system (it de- 
pends only on the observation time t — t w ), and governs 
the dynamics for the shorter observation times. Many re- 
sults in spin glasses showed that the most appropriate form 
for the decay is a power law with a small exponent, a w 
0.1. This behaviour is consistent with the quasi-equilibrium 
noise power spectrum, which varies as l// 1 ^". As this part 
is stationary, it should behave as in the equilibrated system: 
FDT should hold between the stationary part of the relax- 
ation and the corresponding part of the autocorrelation, as 
shown in the section lTTl 

- The second one is non-stationary and decays approxi- 
mately as a stretched exponential of the ratio t/t w . This 
means that if t w — ► 0, this part tends to be instantaneous, 
and if t w — > oo, it becomes infinitely slow. This contribu- 
tion can be rescaled using a re-parametrised (effective) time 
A(i). When plotted versus the effective time difference, all 
the non-stationary contributions measured with different 
waiting times merge very satisfactorily in one curve 1131 . 
showing that the same kind of dynamics persists along the 
whole experimental time range, as in Fig. [6] Here, FDT 
cannot be of any help to link response and stationary parts. 

The limit in the definition of the effective temperature 
(equation |3 means that the two contributions must be well 
separated, i.e., the stationary dynamics must become negligi- 
ble before the ageing one starts to be effective. This situation is 
referred as the "time-scale separation limit", and the evolution 
of any dynamic quantity should present a plateau (in log-scale 
of time) separating the stationary dynamics at short times from 
the ageing one. Experimentally, this clear separation of the two 
contributions is not observed. 

In section|2] a setup allowing the measurement of magnetic 
fluctuations and the response to the conjugate field will be de- 
scribed, and it will be shown that it allows an absolute mea- 
surement of the temperature, following Eq.[2 In section|3] the 
experimental procedure to study the ageing regime of a spin- 
glass using this tool is described. The results allow to check the 
validity of the effective temperature concept, following Eq. |2j 
and analysed according to various models in section|4] 

2 An FDT-based thermometer 

In this section, a new experimental setup, designed to measure 
quantitatively the relations between fluctuation and response in 
magnetic systems will be described. It will be shown that this 
setup works in fact as an absolute thermometer. 

Using FDT, expressed as in Eq. ^ f° r instance, any sys- 
tem allowing a quantitative comparison between thermal spon- 
taneous fluctuations of an observable and the response to its 
conjugate field allows an absolute determination of the temper- 
ature. The new experimental setup developed for the studies 
of spin-glasses is first of all an absolute thermometer, which 
should allow a determination of the thermodynamic temper- 
ature of any equilibrated magnetic system down to very low 




Fig. 1. Schematic of the basic circuit for noise measurement. In order 
to maximise the coupling factor between the sample and the Pick Up 
coil, a long cylinder (4 cm long, 4 mm wide) is used. The third order 
gradiometer being 2.2 cm long 5.5 mm wide, this size insures almost 
the best possible coupling factor, as any contact between the PU-coil 
and the sample must be avoid to allow the temperature regulation. 

temperature. In a setup completely dedicated to low temper- 
ature measurements, the lowest temperature to be measurable 
should be below the milliKelvin range. 



2.1 Noise measurements 

The protocol of spontaneous magnetic fluctuations measure- 
ments is quite simple: a thermalised sample is introduced in 
a pick-up coil (PU), itself part of a superconductive circuit 
involving the input coil of a SQUID-based flux detector (see 
Fig.ID. 

Materially, the sample of cylindrical shape with diameter 
and length 5 mm and 40 mm respectively is contained in a 
cylindrical vacuum jacket, part of a 4 He cryogenic equipment. 
The PU is wound on the jacket. The sample itself is contained 
in a cylinder made of copper coil-foil whose upper part is 
a copper sink with thermometer resistor and heater resistor, 
thermally connected to the 4 He bath by a flexible copper link, 
thus allowing temperature regulation above 4.2 K. The vacuum 
system involves a charcoal pump thermally connected to the 
4 He through a thermal impedance. When cold, it insures good 
thermal insulation; if heated, it allows to inject 4 He exchange 
gas thus thermally connecting the whole sample to the helium 
bath temperature. 

The difficulties of the measurement lie in the extreme 
weakness of the signal of the fluctuations and the strong re- 
sponse to external excitations: the typical amplitude of mag- 
netic fluctuations in our macroscopic CdCr1.7Ino.3S4 sample 
corresponds to the response to a magnetic field about 10~ 7 G. 
Several magnetic shields (^-metal and superconductive) are 
used in order to decrease the residual field at a level of order 
1 mG, and to stabilise it. Furthermore, the PU is built with 
a third order gradiometer geometry, which strongly reduces 
the sensitivity to the time variations of the ambient fields. 
In such conditions, and because of the extreme sensibility of 
SQUID measurements, a satisfactory signal/noise ratio can be 
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easily obtained at short time-scales, corresponding to corre- 
lation measurements with time differences of few seconds. In 
order to study a glassy system in the deep ageing regime, such 
time-scales are not enough: one needs to measure the time 
autocorrelation with time differences up to several thousand 
seconds. To suppress spurious drifts of the measuring chain, 
further precautions are then needed: stabilisation of the helium 
bath to avoid drifts of the SQUIDS sensor, stabilisation of the 
room temperature to avoid drifts in the ambient temperature 
electronics, etc. 

It should be emphasised that the use of the third order gra- 
diometer in this experiment is quite different from the most 
common use. Usually, gradiometers are used in magnetome- 
ters where the sample is small compared to the gradiometer 
size, and is placed in a non-symmetric position; an homoge- 
neous field is established, and the unbalanced flux due to the 
magnetisation of the sample is recorded. 



2.2 Response measurements 

Already several years ago, a comparative study of the mag- 
netic fluctuations and the conventional magnetic response was 
done using a setup similar to the one schematically described in 
figure ^ 1 14 1 . This comparison could not be made quantitative 
with a satisfactory accuracy: the coupling factor of the sam- 
ple to the detection system depends on the PU geometry and 
is quite different in the noise setup and in a classical magne- 
tometer with homogeneous field. In order to be able to compare 
quantitatively the results of both kinds of experiment, one has 
to eliminate the effect of this geometrical factor in the compar- 
ison. This can be done only if the coupling factor is the same 
in both experiments. The way to achieve this can be illustrated 
very simply. The fluctuations of the magnetisation are recorded 
through a PU coil with a given geometry. According to the reci- 
procity theorem, the flux fluctuations induced in the coil are the 
fluctuations of the scalar product of local magnetisation MdV 
by the local field h produced by a unit of current flowing in the 
PU coil: 

<£= / M.hdV 

J sample 

In our setup, the measured fluctuating observable is the flux 
in the PU-coil. The conjugate quantity of this flux is the current 
flowing through the coil, and thus, the magnetic field conju- 
gate of the sample magnetic moment is the field produced by 
the PU-coil itself. If this field is used as exciting field, then 
the fluctuation-dissipation relation should remain the same for 
the macroscopic quantities as for the microscopic ones. This is 
strongly different from the situation where one tries to compare 
the results of noise measurements to the results of classical re- 
sponse measurements done in an homogeneous field: then the 
coupling factor in both measurements has to be evaluated. A 
way to use the PU-coil as field generator is the following. A 
small coil coupled to an excitation winding with mutual induc- 
tance M is inserted in the basic superconductive circuit (see 
Fig.Ea). 




Fig. 2. Schematic of the FDT circuits. TOP : Basic FDT circuit. BOT- 
TOM : The bridge configuration used. 

2.3 Absolute thermometer 

2.3.1 Fluctuation Dissipation relation 

Here we will show that, for any given equilibrated system, the 
validity of FDT on microscopic quantities results in the validity 
of an "effective FDT" on measured quantities. The factor K 
which appears is setup dependent, but sample independent. 

When a magnetic sample is inserted into the PU coil, by the 
reciprocity theorem, a moment m, at position induces in the 
coil a flux 5<P = m;hj . Therefore, the flux in the coil due to 
the sample is given by 

i /i 

where fj, indexes the spin components: fj, = {x,y,z} for 
Heisenberg spins, \i = {z} for Ising ones, etc. We suppose that 
the medium is homogeneous, the components of the fluctua- 
tions are statistically independent and their spatial correlations 
are much smaller than the scale of the PU: 

(m?(l/)m%(t)) = (m(t')m(t)) 6 l3 S^. (4) 

Then, the flux autocorrelation is given by 

(mm) = E E h f m*x*)> = w , *). (5) 

The flux autocorrelation in the PU is thus the averaged one site 
moment autocorrelation per degree of freedom C(t' , t), multi- 
plied by the coupling factor Q determined by the geometries of 
the PU field and of the sample. 
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On the other hand, the impulse response function of one 
moment in the sample is given by 



R^(t',t) = -^ r = R(t',t)S ij 6^, 



(6) 



where R{t! ', t) is the averaged one site response function of the 
sample. If a current i is flowing in the coil, the flux on the coil 
due to the polarisation of the sample reads 



*(*) = EE^ 2 / iZ( * , '* )i( * ,)tft '' 



(7) 



Thus, the response function of the flux due to the sample in the 
PU circuit is 

Rp(t',t) = QR(t',t). (8) 

The same coupling factor Q = J2i 2 determines the val- 

ues of correlation and response of the flux due to the sample. 

Note that the term h is the value of the internal field in the 
sample, due to a unit of current flowing in the coil. Therefore, 
Q corresponds to the same demagnetising field conditions in 
both measurements. Actually, Q is time dependent, since the 
internal field is h = hofj,(t', t) where ho is the field term gen- 
erated by the coil in vacuum and [i(t', t) is the time dependent 
sample permeability, but the important point is that Q(t' ', t) is 
exactly the same in both experiments. 

The above derivation is done in the context of a magnetic 
system, showing that the measured quantities represent those 
used in theoretical work, in which the single-site autocorrela- 
tion and response functions are computed, and averaged over 
the sample. Incidentally, an equivalent derivation could be done 
for any system with magnetic response, for instance the eddy 
currents in a conductor, with the same result: the coupling fac- 
tors are the same in the fluctuations and the response measure- 
ments. 

In the basic measurement circuit, Fig.|2]a, the total flux im- 
pulse response of the circuit to the current i(t') flowing in it 
is 

Rdt', t) = E L8(t - t') + Q{t', t)R{t', t), (9) 

where L is the total self inductance of the circuit. Flux con- 
servation in the (SC) PU circuit leads to 



&exc(t) 



R L (t',t)i(t')dt' = 0, 



(10) 



where ^ ea;c (i) = MI{t) is obtained by injecting a current I{t) 
in the excitation winding. The conjugate variable of the circuit 
current i is the flux <? exc injected by the excitation coil. In the 
case of an ergodic sample, it is easy to show that, once FDT 
applies to the fluctuations and response of the flux induced in 
the PU, it applies also to the fluctuations and response of the 
current flowing in the circuit. Thus, 



k B l 



(11) 



The SQUID gain is G = V s /i. Thus, if a current I(t) = 
Iq(1 — 0(t)) is injected in the excitation coil, the relaxation 



of the SQUID output voltage is related to the autocorrelation 
function of its fluctuations by: 



(12) 



where K = -jj^kb ■ The system is an absolute (primary) ther- 
mometer since, by measuring both the response voltage to an 
excitation current step and the autocorrelation of the voltage 
free fluctuations, it allows a determination of the temperature 
whose precision (once a sample with large signal is chosen) 
depends only on the precision of the determination of the ex- 
perimental parameters Iq, G and M. 

The main drawback of the elementary measuring circuit de- 
picted above is that the response to an excitation step involves 
the instantaneous response of the total self inductance of the 
circuit (first term in the right hand side — R.H.S — of Eq.|9). In 
our case, both the susceptibility of the sample and the coupling 
factor Q are weak. The quantity to be measured, — the second 
term in the R.H.S of Eq. |9] — , represents a few percent of the 
first one. Thus, a bridge configuration as depicted in Fig. [2]b 
has been adopted. Now, the main branch involving the sam- 
ple is balanced by an equivalent one without sample. This sec- 
ond branch is excited oppositely, in such a way that when the 
sample is extracted from the PU, there is no response of the 
SQUID to an excitation step. When the sample is placed into 
the PU, the response of the SQUID is determined only by the 
response of the sample. Nevertheless, now, the loop coupling 
factor of the sample to the SQUID involves different self in- 
ductance terms in both measurements, and one gets 



G L + 2L S 
K = ~r~rz k B 



MI a 



(13) 



where Lq and Ls are the self inductances of the PU and of 
the SQUID input respectively, and the effect of the sample has 
been neglected in the value of Lq- This adds sources of error 
on the calibration since the self inductance values are difficult 
to determine precisely. 



2.3.2 Calibration. 

The circuit as described above is a thermometer, allowing 
the determination of the temperature, ksT. Nevertheless, it 
involves several home-made coils whose self-inductances can- 
not be determined in their experimental environment without 
large errors. This dramatically limits the precision on the de- 
termination of the temperature. A calibration was thus needed. 
For this, the fluctuations and response of a high conductivity 
copper sample were measured in the setup. This high purity 
(99, 999%) sample has a very low residual resistivity at low 
temperature, obtained by annealing at high temperature in oxy- 
gen atmosphere, thus reducing the density of magnetic residual 
impurities. The sample has a cylinder shape, 5 mm wide and 
4 cm long. It was thermalised at the temperature of the boiling 
4 He at normal pressure (4.215 K). 

Since this equilibrated system is stationary, one can use 
standard Fast Fourier Transform algorithms in order to com- 
pute the autocorrelation function from a single record. The 
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Fig. 3. Measured relaxation (a) and autocorrelation (b) function for 
the copper sample at 4.2 K. (c) FD-plot, relaxation versus response, 
the observation time r = t—t w being used as parameter: the observed 
linear behaviour allows calibration of the system as thermometer, this 
slope being proportional to 1/T, independently of the sample. The 
observed deviation from the linear behaviour at the shorter times is 
due to the effect of the low-pass filtering of the excitation — which 
does not affect the correlation measurement. 

average of the obtained autocorrelation function over many 
successive records allows to reduce the noise level. 

The relaxation function is obtained as the response to a field 
step at t w , and is only a function oft — t w . As the system does 
not have remanent magnetisation (the eddy currents vanish in 
a finite time, a few tenth of a second), the limit value of the 
response function is zero. 

The measured relaxation and autocorrelation of SQUID 
voltage are plotted in figure[3ja) and (b) versus the observation 
time. The fluctuation-dissipation diagram (FD-plot) is obtained 
by plotting the relaxation versus autocorrelation, using the ob- 
servation time as parameter, in figure |5Jc) . The observed 
linear behaviour is consistent with the FDT. As previously 
shown, the slope between the relaxation and the response is 
sample-independent, and proportional to 1/ KT. The measure- 
ment on the copper sample, at a well known temperature allows 
thus to determine the factor K. From it, we can determine the 
temperature of any sample placed in the gradiometer from the 
value of the slope of the measured relaxation versus correla- 
tion curve. This can be applied to any equilibrated system. For 
glassy systems, it should allow an experimental determination 
of the effective temperature. 

3 Fluctuation-dissipation relations in a 
spin-glass sample 

The aim of this work is the study of the fluctuation-dissipation 
relation in a spin glass. In this section, we will emphasise the 
peculiarities of this measurement, and then describe the proce- 
dure used to make the analysis quantitative as well as the limits 
of this procedure. 



3.1 Sample 

The knowledge acquired in previous magnetic noise investiga- 
tions on spin glasses was very helpful to choose a good candi- 
date for the present study. First, eddy currents in metallic sam- 
ples produce noise, as in the copper sample used for calibration. 
This noise can be measured, but not directly related to the spin 
dynamics. In order to avoid this drawback, an insulating spin 
glass sample was chosen. 

Measurements on CsNiFeFg have shown that this com- 
pound has a stronger signal, and then a better signal to noise 
ratio than any other insulating spin glass 1151 . Nevertheless it 
has a very strong ferromagnetic value of the average interac- 
tion and its behaviour is far from the "standard" spin glass be- 
haviour. 

CdCr2_2xIri2xS4 was also extensively studied experimen- 
tally, by classical susceptibility, magnetic noise, neutron scat- 
tering II 1 61 1 71 1 81 1 91 . In this series of compounds, the magnetic 
ions are Cr 3+ , with low anisotropy. The coupling is ferromag- 
netic between first neighbours, and anti-ferromagnetic between 
the second ones. The substitution of Cr 3+ by the non-magnetic 
In 3+ increases the relative importance of anti-ferromagnetic 
coupling as compared to the ferromagnetic one. The random 
dilution introduces disorder and frustration, the basic ingre- 
dients leading to spin-glasses. In the studies of the spin glass 
state, CdCr1.7Ino.3S4 is the preferred compound in this fam- 
ily. The high concentration of Cr 3+ allows this sample to have 
a strong signal, but it is not high enough to reach the percola- 
tion of the ferromagnetic order. With decreasing temperature, 
finite sized ferromagnetic cluster formation is observed. Close 
to T g , these clusters are rigid, and the interaction between 
them is random, with a weak anti-ferromagnetic average. This 
clustering greatly increases the signal, as the noise power of 
N = Nq/ti ferromagnetic clusters of n spins is y/n times 
stronger than the one of Nq individuals spins. 



3.2 Experimental details 

Glassy systems are not stationary; their dynamics depends on 
two times, both referred to a crucial event, the "birth" of the 
system. In spin glasses, the birth time is best defined by the 
time at which the final temperature is reached, as soon as the 
end of the cooling procedure is fast enough 1201 . A cooling pro- 
cedure based only on driving the sample holder sink tempera- 
ture would introduces non-negligible temperature gradients if 
the cooling or heating rate is too high. This would lead to a 
distribution of ages over the sample. In order to obtain a more 
homogeneous temperature, the cooling procedure is as follows: 

- first, the temperature is slowly decreased from a reference 
temperature T re f above T g to a temperature T\ w T m + 
3 K, where T m is the working temperature; 

- then, by heating the charcoal pump, a small amount of He 
gas is introduced, allowing a quick and homogeneous cool- 
ing; 

- finally, the charcoal pump heating is switched off and the 
vacuum surrounding the sample is restored, allowing the 
temperature regulation at T m . 
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Fig. 4. A typical thermal history of the sample for a 4, 500 srecord at 
10 K. In inset, detail on the crucial part, the last 3 K's cooling. 



The first step (slow cooling, approx. 0.25K.min _1 ) does not 
introduce severe temperature gradients, as the cooling is slow 
enough. This step cannot be avoided, as the fast cooling by 
exchange gas can only be used to decrease the sample temper- 
ature by few Kelvin without introducing too strong perturba- 
tions in the helium bath. Previous studies on CdCr2-2xIn2 X S4 
show that the dynamics is governed by the second, fast, cooling 
step, at least on timescales shorter than 20 , 000 s 1 20 1 . This sec- 
ond step is obtained by heating the charcoal pump during few 
seconds. As the gas surrounds the sample, the resulting cool- 
ing is homogeneous. When the heating is stopped, the charcoal 
absorbs the gas back, allowing the regulation of temperature. 
The amount of gas used and the duration of heating are ad- 
justed in order to lower the temperature exactly down to the 
working temperature. This allows to cool the sample by 3 K in 
less than 30 s, keeping the temperature gradient negligible. The 
birth time is taken as the instant when the sample temperature 
reaches T m + 15 mK. This allows a precise, and reproducible 
determination of it. 

The measurement of the relaxation is straightforward. A 
DC current is applied to the excitation coils at high tempera- 
ture, before the beginning of the quench procedure, and switched 
off at t = t w . The relaxation is then recorded for t > t w . The 
signal is measured before applying the excitation: this deter- 
mines the zero baseline of the measurement. After relaxation, 
the sample is re-heated to the start temperature in order to 
check the stability of the baseline. Thus, in the measurement, 
both the zero and field cooled (FC) levels are known. 

Recording the fluctuations is even simpler, at least in princi- 
ple: no field is applied, the spontaneous fluctuations of the sig- 
nal are just recorded from the end of the quench procedure and 
during a long enough time to be able to compute all the desired 
C(t w , t). However, as the system is not ergodic in the ageing 
regime, the autocorrelation of the signal cannot be evaluated 
from a single record, as for an equilibrated sample. In order 
to compute the autocorrelation function, an ensemble average 
has to be done on successive equivalent records, each one initi- 
ated by a quench. This does not only increase dramatically the 



length of the experiment, but also the difficulties of the acqui- 
sition, the ideal acquisition conditions having to be kept during 
months instead of hours. This complication has nevertheless 
an advantage: by averaging, it allows to make the separation 
between a systematic spurious signal and the signal of the fluc- 
tuations. In our results, the systematic signal is of the same or- 
der of magnitude as the fluctuations themselves. It corresponds 
to the drift of the SQUID due to the continuous decrease of 
the He-level, and to the global response of the sample to the 
residual field during the cooling procedure. The average of the 
signal over records gives the zero, and the sample fluctuations 
signal is then given by: 

m(t) = M(t) - (M(i)) (14) 

The autocorrelation is then evaluated from its definition: 

d(t,t w ) = (m(t w )m(t)) (15) 

In order to obtain a small statistical error, a very huge number 
of records is needed. As each record length is about few hours, 
the number of records is limited to about 300, and the aver- 
age over the records is not enough to obtain a satisfactory ratio 
between the statistical error and the signal. As the autocorre- 
lation function should evolves smoothly for both variables, t w 
and t, the autocorrelation function computed following eq. ^] 
is averaged over small time intervals of both variables : 



C'avg(to,t Wo ) — Cl(*,iw)|t„e[t TOO ±e 4lO0 ],te[t o ±e t0 ] (16) 

with 

e t = 0.05i<t. (17) 

The criterion used (Ea. ll7t is a compromise between the need 
of statistics in order to obtain a low enough statistical noise, 
and the requirement to keep e as small as possible to be able to 
capture as precisely as possible the non-equilibrium dynamics. 

3.3 Correlation offset. 

In principle, our experimental procedure, involving many reali- 
sations of the same experiment, allows to compute the autocor- 
relation function of the magnetisation following its exact defi- 
nition, and thus exactly. Nevertheless, in reality, this would be 
the case only if external sources of noise were negligible, not 
only in the correlation time-scale under study, but also in the 
time-scale of one complete record. This means that the external 
noise should be controlled not down to 1 mHz as in our experi- 
ment, but at least down to frequencies as low as few 0.01 mHz, 
which is quite impossible. The result is that the computed cor- 
relation C(t w , t—t w ) contains an offset practically independent 
on t — t w but randomly dependent on t w . 

As the setup is a calibrated thermometer, the temperature 
can be extracted from the derivative of the x(C) curves in 
experimental units. In order to obtain the FD-plot, this is not 
enough. For a normalisation of the data, one needs to know 
the zero reference level of the response and of the correla- 
tion. In the case of the copper sample data, where the eddy 
currents producing the signal have a finite and experimentally 
accessible lifetime, this calibration is trivial: relaxation and 
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autocorrelation functions decrease to zero after a few seconds. 
In the spin glass case, normalisation of the relaxation is simple 
since the zero level and the FC level are determined during the 
measurement. For the correlation, things are not so easy. 

As seen above, at a given temperature, the correlation 
curves for different t w are shifted between each other by a ran- 
dom unknown offset. Nevertheless, it is possible to normalise 
the data by taking as the origin of each curve the square of 
the measured value of the first point, C avg (t w ,t w ). Due to the 
elementary measurement time constant, this term corresponds 
to an average over t—t' about 10 -2 s, i.e., a range of (t — t') /t' 
corresponding to the stationary regime where all curves must 
merge. Thus, the following quantity is computed: 

The "individual" offset is now replaced by a "global" one, Co, 
which should apply simultaneously to any measurement done 
at a given temperature. 

Then, the best way to normalise our data could be to ex- 
tract C(t,t) from some other measurement, and to be able 
to convert it in the "experimental units". Neutron diffraction 
experiments are now under way in order to measure this quan- 
tity. C(t, t) can also be extrapolated from high temperature 
measurements (above T g ) to low temperatures (below T g ), or 
deduced from some other quantities. Then a complete — but 
model-dependent — determination of the autocorrelation can 
be obtained, allowing to obtain the FD-plot. Anyway, even if 
the hypothesis used to obtain this complete determination of 
the autocorrelation were not realistic, some characteristics of 
the FD-plot would not be affected. The temperatures, effective 
or not, are measured from the slope between relaxation and 
correlation in the experimental units. They will not be mod- 
ified by the normalisation procedure, whose effect is just to 
suppress an offset. 



4 Discussion 

In this section, the results of the measurements done at several 
temperatures in CdCr1.7Ino.3S4 will be analysed following 
the line of the method described above. 



4.1 Raw measurement 

Figure|5]displays the values of x( t w,t)/xFC = l-<r(t w ,t)/XFC 
plotted versus C(t w , t) for several values of t w and using t — t w 
as parameter. The three graphs correspond to the three temper- 
atures 10, 13.3 and 15 K. A first observation is that a linear 
regime exists between relaxation and correlation for all the 
temperatures and waiting-time investigated. This regime cor- 
responds to the shorter observation times. In the figures, the 
straight lines represent the FDT slope as calculated from the 
values of calibration factor K and of the temperature: in this 
regime, the relation between relaxation and correlation fol- 
lows the FDT. Thus, this regime can be extrapolated from 
the shorter experimental observation-time down to the mi- 
croscopic time-scale. This extrapolation at short time should 
reach the starting point of the FD-plot: C(t w , t w ) corresponds 



t° x(tw, tw) = 0. As no long term memory is observed in spin- 
glasses, C(t w , +00) = should also correspond to x{tw, +00) = 
Xfc, but the extrapolation to this point is not obvious at all, as 
the (unknown) ageing regime should be extrapolated. 

4.2 Scaling procedure 

The raw results show a waiting time dependence which can be 
easily explained. The main theoretical predictions correspond 
to the approach of the limit t w — > 00, C(t w , t) = C (WEB). 
In this case, the stationary and the ageing parts of the dynam- 
ics evolve on distinct time-scales, yielding a separation of both 
dynamics. Experimentally, such a separation is not accessible 
since the waiting times are finite. In order to separate both part 
of the dynamics, the entanglement between both part should 
be described. The simplest way to combine these two contri- 
butions is to add them. For instance, if one considers the relax- 
ation (7 observed after a unitary field step at time t w : 

cr(t w ,t) = (l-A).a stat (t-t w ) + 

A.a ageing (X(t) - X(t w )) (18) 

In this equation, all the different a are normalised to unity. 
This relation is obviously valid in the limit of separation of 
time-scales, and is the most commonly used in theoretical ap- 
proaches, but it is counter-intuitive as shown by the following 
two thought-experiments: 

- In the first one, a glassy system is quenched at a temper- 
ature below T g at time t = 0, and at t = + , a field is 
applied. This experiment is usually though as being equiv- 
alent to the "Field-cooled" procedure, in which the field 
is applied before the quench. Experimentally, the Field 
Cooled magnetisation is strikingly stable. In the additive 
formulation, the predicted behaviour is the following: first, 
an instantaneous variation due to the ageing part, and then a 
slow variation up to the FC value, as the system approaches 
equilibrium. Thus, the field-cooled magnetisation should 
not be as stable in time, as what is observed experimentally. 

- The second problem arises when thinking about some fi- 
nite t w experiments, but with (very) huge time differences, 
t — ► 00. Ageing and stationary parts are usually described 
as stretched exponential (with characteristic time of order 
t w ) and power-law with small exponent respectively. For fi- 
nite t w , after a finite time, the only remaining contribution 
to the dynamics would come from the stationary part, and 
FDT would be recovered. 

If the time-scales are not well separated, it seems intu- 
itively that the two different contributions must be more en- 
tangled than the result of a simple addition. Another (more 
realistic) way to mix the two parts together is to consider a 
"multiplicative" combination, which can be written as: 

a{t w , t) = [(1 - A).<r stat (t - t w ) +A]x 

O ageing (A(t) - X(t w )) (19) 

In the limit of time-scales separation, equations 1181 and[T9l 
are equivalent. Moreover, it is easy to show that the problems 
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C(t w ,t)-Co [V 2 pp ] C(t,t w )-C [V 2 pp ] 




0.0 0.1 0.2 0.3 0,0 0,2 0,4 0,6 0,8 1,0 



C(t,t w )/C(t,t) C(t,t w )/C(t,t) 

Fig. 5. Raw results (full symbols) and ageing part (open symbols) deduced from the scaling analysis for the three investigated temperature, 
T — 10, 13.3, 15 K. The straight lines have the FDT-slope obtained from the copper calibration for each temperature of the formalisation bath, 
and start all from the point (u/xfc = 1, C/C(t, t) = 1). The correlation offset is adjusted following the "PaT hypothesis" (cf section l4,3,3> 
The different curves span the waiting times studied: o : t w = 100 s, A : t w = 200 s, V : t w = 500 s, O : t w = 1000 s, < : t w = 2000 s 
and, only for T = 13.3 K, > : t w = 5000 s, : t w = 10000 s.The last plot represent the complete FD-plot for each previous measurement 
measurements reported on the same graph, the smooth curve corresponds to equation l26l with an exponent B — 0.47. 



raised in the two cases discussed above disappear. Thus, as this 
formulation is at least as justified as the "additive" one, and 
seems less counter-intuitive to the authors, it will be preferred 
in the following. 

Experimentally, separation of timescales is not accessible 
since the waiting times are finite. In order to separate both part 
of the dynamics, a scaling analysis, as described by equation 
1 191 and illustrated by Fig. [6] is applied on both relaxation and 
correlation measurements, within the following constraints: 

(i) In the relaxation, the stationary part is described by a 
power-law decay; its exponent a is extracted from the 
decay of the noise power-spectra measured on the same 
sample, at the same temperature, in the quasi-equilibrium 
regime obtained after a very long waiting at the working 
temperature (typically 15 days) 1121 1 . 



(ii) In the non-stationary regime, the effective time is given by 
A = j—. The value of the sub-ageing exponent, i.e., [i is 
in the range 0.85 — 0.9 III 61 . /x and the relative amplitude 
of the ageing part, i.e., A, are chosen to obtain the best 
rescaling of the relaxation curves once the stationary part 
has been subtracted. 

(iii) The stationary part of the autocorrelation function is then 
evaluated using FDT and the stationary part of the relax- 
ation. The ageing part of the autocorrelation can then be 
deduced from eauation[l"8lorl 1 91 

The resulting FD-plots are displayed in open symbols in the 
diagrams of figure [5] for the ageing part (by construction, the 
stationary part follows the FDT line). For all the investigated 
temperatures, the ageing part starts with a slope very close to 
the FDT one. This may reflect the imperfections of our decom- 
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0.1 1 10 100 1000 
0.4 i 




Fig. 6. Relaxation (Top) and autocorrelation (BOTTOM) functions 
recorded at 13.3 K. The different curves correspond, from bottom to 
top, to t w = 100, 200, 500, 1 000, 2 000, 5 000 and 10 000 s. In insets, 
the respective ageing parts, deduced by the scaling analysis developed 
in section l4~2l 



position between the stationary and the ageing parts. Using the 
additive scaling, this effect is even more pronounced. 

The multiplicative scaling has a main drawback: it can 
not be used without the knowledge of the amplitude of the 
concerned physical quantity. As previously discussed in sec- 
tion 13.31 for the correlation this value must be determined 
indirectly, and is model dependent; thus the scaling is also 
model-dependent. The additive form of the scaling can be 
applied without any amplitude parameter. However, the de- 
pendence of the scaling on this parameter is weak, the results 
obtained from different models are indistinguishable from each 
other. 

The ageing part of the correlation function is found to fol- 
low remarkably the scaling-law used for the ageing part of the 
response, resulting in a very weak systematic t w -dependence of 
the FD-plots. Thus, the equation[T8l(or[T9l can be written for 
the correlation, with the amplitude parameter A replaced by 
qEA, the usual Edwards- Anderson order parameter 1221 . which 
is defined as the remaining part of the autocorrelation for an 
equilibrated spin glass after an infinite waiting time. As a con- 
sequence, the FD-plots are determined by a single parameter, 
the effective time difference: the limit FD-plots, corresponding 
to the ideal separation of the stationary and ageing regime is in- 
dependent of the age of the system, as supposed in theoretical 
works. 



4.3 Comparisons with some models predictions 

Depending on the models, some remarkable features of the FD- 
diagrams are predicted. The analysis of the FD-diagrams may 
help to check the validity of the models used to interpret the 
glassy behaviour found in CdCri. 7I110.3S4. 

4.3.1 Domain growth 

In domain growth models, as in any replica symmetric mod- 
els, the FD-plot should be quite simple in the limit of time- 
scale separation. For an infinite waiting time t w , the quasi- 
equilibrium relaxation should go down to zero. A FDT-behaviour 
should then describe all the response, governed by the single- 
domain response and the scaling approach used in this paper 
should give (qEA 7^ 0, A = 0). Thus, the remaining part of 
the diagram, an horizontal line, should correspond to an infi- 
nite effective temperature 1238241. This description does not 
coincide with the FD-diagram shown in figure|5] even after the 
separation of time-scales obtained by scaling. 

Anyway, a more refined approach as in 1251 is not excluded, 
in which dynamics is described introducing a crossover region 
in between the quasi-equilibrium region (C > qEA) and a 
purely ageing region characterised "dynamical order parame- 
ter" qu for C < qo < qEA- This approach could explain the 
"early" departure from the FDT regime. This departure should 
be t w -dependent, but this dependence may be hidden by a too 
small range of waiting times explored (as well as a too weak 
exploration of the ageing regime). 

4.3.2 1-step replica symmetry breaking 

In CdCr1.7Ino.3S4, one of the best realisation of an Heisen- 
berg spin-glass, it has been shown that the scenario of the chi- 
ral spin glass could be relevant 1261271 . Such model belongs to 
the 1-step replica symmetry breaking (1-RSB) models family 
H28I29I30I . 

In 1-RSB case, the ageing regime is described by a unique 
effective temperature, finite and strictly greater than the ther- 
malisation temperature. By considering a normalised FD-plot, 
it is easy to show that the value of qEA can be deduced from the 
value of T, T e f / and 7 = the ratio between the stationary 
and the ageing part of the relaxation, which is experimentally 
accessible: 



As the experimental setup is a calibrated thermometer, it 
allows an absolute determination of the temperatures. The de- 
termination of T and T e f f extracted from the slope of the sta- 
tionary and the ageing part respectively allows the complete de- 
termination of the offset C . The obtained values of T e ff and 
qEA are reported in tabled The separation by scaling between 
stationary and ageing part being far from perfect, the uncer- 
tainty on the determination of T e ff, and consequently on qEA 
is quite large. The choice of a scaling procedure also influences 
the results (the previously reported value for T e f f sw 30 K for 
a thermalisation temperature of 13.3 K was obtained by an ad- 
ditive scaling analysis of the data 11311 ). The results of table 
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T[K] 


Teff[K) 


1 — step 

Iea 


PaT 

Iea 


10 


28 ±6 


0.65 ± 0.05 


0.63 


13.3 


50 ± 10 


0.45 ± 0.05 


0.37 


15 


80 ±20 


0.36 ± 0.07 


0.21 



Table 1. Values of T e // as obtained in an 1-step replica symmetry 
breaking scenario and the corresponding Edwards- Anderson order pa- 
rameter ev (cf section I4.3.2L The values of qEA T 316 deduced 
from the PaT ansatz (cf. sectionl4.3.3>. 



^can be compared with the results of simulations done on a 
weakly anisotropic spin glass model by Kawamura 1321 . In this 
work, it was found that \ depends linearly on C in the ageing 
regime. The effective temperature was found to be indepen- 
dent of the temperature of the thermalisation bath. This inde- 
pendence does not appear in our data, but, maybe, it can be 
due to the extremely low anisotropy used in the simulations. 
CdCr1.7Ino.3S4 is known to be an Heisenberg spin glass with a 
non negligible anisotropy, which as been found to be five times 
stronger than in the canonical AgMn spin glass. 



4.3.3 Continuous replica symmetry breaking 

In continuous replica symmetry breaking (oo-RSB) models, 
as the Sherrington-Kirkpatrick (SK) model 1331341 . the Parisi 
order parameter is a continuous function between and qEA 
1351361 . Links between statics and dynamics imply that the 
corresponding effective temperature is a smooth and not trivial 
function of the autocorrelation : 



T 



eff 



1 

T 



P(q)dq 



(21) 



Then there is no trivial relationship between qEA and the 
measured quantities. It is nevertheless possible to to progress 
further if the studied compound is a canonical spin-glass. In 
these systems, the FC susceptibility is purely paramagnetic 
at high temperature, following equation [22] and below T g , its 
value is temperature independent: 



Xfc(T > T g ) 

Xfc(T < Tg) 



= C(t,t)/k B T 

C(t, t)/k B Tg 



(22) 
(23) 



The lower the concentration of magnetic ions in the canon- 
ical sample, the smaller the probability of spins clustering 
and the better the validity of the above relations 1371 . Thus, 
the value of C(t,t) can be straightforwardly derived from 
susceptibility data. The canonical behaviour is also observed 
or imposed in theoretical models, and known as the Parisi- 
Toulouse Hypothesis 1381 . This "PaT" hypothesis implies that 
the FC response is temperature independent, as observed in 
diluted spin glasses. 

In samples with high concentration in magnetic sites, de- 
viations from the simple canonical behaviour are observed, 
as well as the formation of clusters of spins. At low tempera- 
tures, the response of single spins is no more observed, but the 
response of some rigidly coupled groups of spins. For macro- 
scopic quantities, this is equivalent to the response of fewer 



renormalised spins. In CdCr1.7Ino.3S4, which has a mean cou- 
pling constant strongly ferromagnetic (0 = 100 K 1171391 1. 
this clustering may explain that at low temperature, but above 
T g , the compound behaves as a compound with antiferro- 
magnetic average of couplings. A standard Curie-Weiss law 
description around 30 K gives a mean coupling characterised 
by r* — 9K 1401 . Such a description, with a non-trivial 
0{T) variation, should be associated with a non-trivial, but 
still smooth, function C(t, t; T). As informations on the vari- 
ations of C(t,t;T) are lacking we propose to consider that 
relations l22l and l23l are still valid in the general case. This is a 
strong hypothesis since it amounts to consider that the temper- 
ature variation of xfc is due only to the temperature variation 
of C(t,t;T). One can write: 



C(t,t;T) (T<T) T v 
- 1 -^^ XFC (T<T g )- 



(24) 



A smooth behaviour of C(t, t; T)/UbT around T g can be 
obtained using T* = 17.2 K in the formula . 

This ansatz gives an access to the unknown offset of the 
autocorrelation: i) the starting point of the FD-plot, corre- 
sponding to [C(t w ,t w ;T); a(t w ,t w )] is completely defined, 
ii) C(t w , 00; T) corresponds to the point where the FDT line 

. T* 



reaches the level given by xfc(T < T, 



g) t 



Then the FD- 



graph can be plotted in reduced units as displayed in figure|5jl. 
In this plot, the starting point (C = 1) and the end point 
(C = 0) are temperature independent. Furthermore, it has been 
shown (for some mean field models, and approximately for the 
SK model) that not only these points but also all the ageing part 
of the plot is temperature independent 1411 . It is conjectured 
that this can be still valid in short range models I42f43|44l . 
This property is particularly interesting: it allows, by measure- 
ments at several temperatures, to obtain the whole "master" 
curve describing the ageing behaviour, even if each set of data 
spans a limited portion of the correlation. This feature has been 
already used to obtain the master curve from response data, 
assuming that the separation of time-scales is reached in usual 
susceptibility measurements 145 1461 . 

In the SK model at small C, it can be shown that the master 
curve should behave as 1431 : 



x(C) = VT^c 



(25) 



For correlations close to zero, the slope of the FD-plot, 

X(C) = J C P(q)dq, should asymptotically reach zero, as 
P(0) is known to have a finite value in the continuous RSB 
case. 

Equation |25]can be generalised by allowing any exponent 
different from 0.5: 



X (c) = (i-cy 



(26) 



Such a curve can be easily superimposed to our experimen- 
tal results. Using a coefficient B = 0.47, a single curve can 
describe the ageing regime at all the investigated temperatures 
and for the data close to qea(T). Both the value of qEA(T) 
and the effective temperature close to qEA seem to be well de- 
scribed by Eg. 1261 
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For C <C qEA, at each temperature, the experimental points 
deviate from relation[26] This cannot be due only to the smaller 
signal to noise ratio at the longest timescales, since the effect 
seems to be more pronounced at the highest temperature, where 
the sample signal is the strongest. 

A possible explanation is that the scenario with continuous 
replica symmetry breaking should be associated with a contin- 
uous distribution of timescales describing the system. As the 
limit of separation of timescales is not reached in our results, 
the ageing regime itself is a combination of many timescales. 
The scaling procedure allows to extract the stationary part, 
but not to reach the limit where a full time-scale separation is 
achieved. 

A way to reach the limit could be to iterate the scaling pro- 
cedure on the ageing data to separate the "ageing timescales". 
The ageing regime can be considered as a pseudo-FDT one, 
with a temperature equal to T e ff(qEA)- Such a work on the 
available data is however hopeless, as the separation between 
stationary and ageing part seems obviously already far from 
perfect. 

A scaling can be deduced from equation 1261 143 1, which, 
using <P = Yz~g> can be written as, : 

w _ f A[(l - C)T-*] B for C < q EA {T) , 
X _ \(l-C)T- # iorC>q EA (T). y£l) 

If a power-law can describe the ageing dynamics, then all 
the scaled data should merge along a single line. The best result 
is obtained for B = 0.5, but the cloud of points remains very 
broad, and is not well described by the predicted straight-line 

in the log(T 1 -' i> x) vs log (T~*(l - C)) diagram. 
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10 K, t =200s 
10 K, t w =500s 
13,3 K, t^lOOs 
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T*(1-C) 

Fig. 7. Scaling of the ageing parts of the FDT-diagram following equa- 
tion^] The straight line shows the result predicted by this equation. 



value of several elements of the setup, the setup was calibrated 
by measuring the magnetic fluctuations and response of a high 
conductivity copper sample thermalised at helium tempera- 
ture. This calibrated thermometer was used to determine the 
out-of-equilibrium properties of a spin-glass below the glass 
transition. The autocorrelation function of the spontaneous 
magnetic fluctuations of a well characterised insulating spin- 
glass was investigated in the deeply non-stationary regime. 
Its waiting time dependence can be described by using the 
same scaling as for the response function. The FD-plots clearly 
confirm that the stationary dynamics observed at the shorter 
timescales can be considered as a quasi-equilibrium once, as 
the fluctuation-dissipation relation between autocorrelation 
and relaxation obeys the fluctuation-dissipation theorem. The 
results show clearly that the asymptotic regime, with full sep- 
aration of timescale is not reached. Certain Hypotheses on the 
dynamics are made in order to compare the results with model 
predictions. The deduced scaling analysis allows to extrapolate 
the experimental results to the limit used in theoretical studies 
of weak-ergodicity breaking models. 

The experimental results obtained on CdCri ,7^0.384 dif- 
fer qualitatively from the predictions of any domain-growth 
model. 

The experimental data allow interpretations that are rather 
consistent with predictions from the two replica symmetry 
breaking models under study. As long as the autocorrelation 
cannot be determined completely, both models can be relevant, 
giving slightly different results concerning the value of q E A- 
An independent determination of the characteristic magnetic 
moment of the clusters as a function of temperature is needed 
in order to resolve this ambiguity . 

The possibility of analysing the experimental results on 
the basis of 1-step replica symmetry breaking confirms that, 
at first sight, the chiral model developed by Kawamura could 
be the more relevant one for the Cd&i 7iiio.3S4Compound 
with low anisotropy, supporting the conclusion from D. Petit 
and I. Campbell on this sample 1261 . However, the scatter of 
the data is such that one cannot reject definitely an interpre- 
tation inspired by the mean-field models developed for Ising 
spin-glasses. 
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